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Abstract 

In three-dimensional QED with a Chern-Simons term we study the phase 
structure associated with chiral-symmetry breaking in the framework of the 
Schwinger-Dyson equation. We give detailed analyses on the analytical and 
numerical solutions for the Schwinger-Dyson equation of the fermion propa- 
gator, where the nonlocal gauge-fixing procedure is adopted to avoid wave- 
function renormalization for the fermion. In the absence of the Chern-Simons 
term, there exists a finite critical number of four-component fermion flavors, 
at which a continuous (infinite-order) chiral phase transition takes place and 
below which the chiral symmetry is spontaneously broken. In the presence 
of the Chern-Simons term, we find that the spontaneous chiral-symmetry- 
breaking transition continues to exist, but the type of phase transition turns 
into a discontinuous first-order transition. A simple stability argument is given 
based on the effective potential, whose stationary point gives the solution of 
the Schwinger-Dyson equation. 

ll.30.Qc, ll.30.Rd, ll.lO.Kk, 11.15.Tk 
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I. INTRODUCTION 



This paper is the detailed exposition of our previous letter |^ on chiral-symmetry break- 
ing in (2+l)-dimensional QED with flavors of four-component Dirac fermions (here- 
after called QED3). Such a model of QED3 is chirally symmetric in the absence of a bare 
fermion mass term, mQipilj, in sharp contrast to the (2+l)-dimensional gauge theory with 
two-component fermions, where we cannot define chiral symmetry [@-^. Similarly to the 
four-dimensional case |p, the chiral symmetry of QED3 may be broken spontaneously due 
to the dynamical generation of a fermion mass. However the pattern of dynamical symmetry 
breaking is shown to be qualitatively different from the four-dimensional counterparts. 

Recently, QED3 has found a vast region of application in condensed matter physics as an 
effective theory in the long wavelength (or low-energy) limit of a more realistic microscopic 
model Especially, since the discovery of high-Tc superconductivity and fractional 

quantum Hall effect, a peculiarity of three-dimensional gauge theory, i.e., the existence of a 
Chern-Simons (CS) term 

^e'^-^'^A^d^A, (1) 

has been a basic ingredient in these applications. The CS term gives the gauge field a 
mass without destroying the gauge invariance. And, even if there is no bare CS term 
in the original Lagrangian, a CS term can be generated by radiative corrections in three 
dimensions. In the presence of a CS term, parity is broken explicitly, but the Lagrangian 
still keeps the chiral symmetry when rriQipilj = 0. The parity violation is a very important 
factor for the anyonic model to be a candidate theory of high-Tc superconductivity 
A field-theoretic realization of such an anyonic model consists of fermions interacting with 
an abelian statistical gauge field whose dynamics is governed by a CS term. In contrast to 
the chiral symmetry, it is known that the dynamical breakdown of parity does not occur in 

QED3 fig-g. 

In this paper we pay attention to the breaking of chiral symmetry in the presence of a 
CS term within the framework of the Schwinger-Dyson (SD) equation In momentum 
space, the SD equation for the full fermion propagator S{p) is written as 

S'^ip) = P/.7m -mo-e^ J -^-^-f^S{q)T^{q,p)D^^{k) , (2) 

which should be solved self-consistently and simultaneously with other SD equations for the 
full photon propagator D^^{k) and the full vertex function V^{q,p). In order to solve these 
sets of SD equations actually, however, we adopt some kind of approximations as exemplified 
below. 

There are several quantities we must specify to write down the closed SD equation for 
the fermion propagator S{p), namely, the full photon propagator D^iy{k) and the vertex 
function r^(p, q), and we also have to choose a suitable gauge. In general, both the photon 
propagator and the vertex satisfy their own SD equation, and both depend on the gauge- 
fixing parameter a, if the covariant gauge-fixing Lagrangian ^{d^A^Y is taken. There are 
also additional constraints on these n-point functions, coming from relations such as the 
Ward-Takahashi (WT) identity 
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{p,-q,)T^ip,q) = S-\p)-S-\q). (3) 

Now one of the problems in QED3 is to find a consistent truncation sclieme for the SD equa- 
tion: as a consequence of the WT identity the wave-function renormahzation should be equal 
to the vertex renormahzation, which puts a constraint on the bare vertex approximation. 
By using 

S{p) = [^,p,Aip)-B{p)]-\ (4) 

the SD equation for S{p) is decomposed into two integral equations for A{p), the wave- 
function renormahzation, and B{p). Although with a bare vertex it is not possible to satisfy 
the WT identity exactly, we should make sure to satisfy the constraint A{p) ~ 1 when using 
the bare vertex approximation. 

First of all, we restrict our attention to QED without CS term. In the quenched case, 
in which the vacuum polarization to the photon propagator is completely neglected, the 
simplest treatment is to choose the Landau gauge a = and to take the bare vertex 
r^(p, g) = 7^. This quenched approximation is equivalent to taking the N ^ limit in 
the SD equation for the fermion propagator, since this limit eliminates the radiative cor- 
rection to the photon propagator coming from the internal fermion loops (with species). 
Then the above SD equation leads to no wave-function renormahzation: A{p) = 1, irrespec- 



tive of B{p), see e.g. [^. In this case, it has been shown that quenched QED3 resides 
in only one phase where chiral symmetry is spontaneously broken, in agreement with the 
Monte Carlo simulation of non-compact lattice QED3 pOl . 



In the quenched case, the bare vertex approximation is justified in the Landau gauge 
a = 0, because this choice is approximately consistent with the WT identity and satisfies 
the requirement that the wave-function renormahzation and the vertex renormahzation are 
equal. In gauges other than the Landau gauge however, we can not take the bare vertex, 
since A{p) deviates from one. In order to study the SD equation in a general gauge, we 
must modify the vertex in such a way that the result is gauge-covariant, i.e., independent 
of the gauge-fixing parameter a. This problem persists in the three-dimensional case pT|J22| 
as well as in four dimensions [p^-p5[|. Such a modified vertex has to satisfy the WT identity 



exactly. 

Now we can raise the question whether the chiral symmetry restores at a certain non- 
zero value of A^. In the presence of the one-loop vacuum polarization (A^ ^ 0) the simple 
procedure mentioned above can not be applied, since there is no simple gauge choice a such 
that A{p) = 1. Nevertheless, under the bare vertex approximation in the Landau gauge, 
Appelquist, Nash and Wijewardhana (ANW) have shown that there is a finite critical 



number of flavors Nc, above which the chiral symmetry restores. Based on a leading-order 
1/A^ expansion (e^ = 8a/N with a being kept fixed) for the vertex and the photon 
propagator in the SD equation, they gave the critical value Nc = 32/7r^ = 3.2. Furthermore, 



Nash p8[ claimed that the leading-order gauge-invariant critical number of flavors is given 
by Nc = 128/37r^ and that, when 1/N^ corrections are included, Nc coincides approximately 
with this result. 

However, such a simple treatment of the vertex function and the wave-function renor- 
mahzation function A was criticised by Pennington and Webb and Atkinson, Johnson 
and Pennington |30|. They claimed that, if the vertex is correctly improved, using the 
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WT identity, and the coupled equations for A{p) and B{p) are solved together, the finite 
critical number of fiavors does not exist. This implies that the chiral symmetry is sponta- 
neously broken in QED3 for all values of A^, i.e., Nc = oo. On the other hand, the Monte 
Carlo simulation of lattice non-compact QED3 seems to support a finite critical fiavors: 
Q Ai"c ~ 3.5 ± 0.5, see also However, it is numerically very difficult to confirm the 

exponential decrease of the dynamical mass for increasing found in ||29|j3^ . 

The origin of this controversy stems from the fact that in the unquenched case A{p) = 1 
can not be deduced as a simple consequence of the SD equation by choosing the Landau 
gauge, in sharp contrast with the quenched case. Nakatani has proposed to use a nonlocal 
and momentum- dependent gauge function, instead of the usual gauge-fixing term, in order 
to keep A{p) = 1 and thus to overcome the inconsistency one has using the bare vertex 
approximation in the Landau gauge. Actually the nonlocal gauge found in p3| , |3^ , p!6| can play 
exactly the same role as the Landau gauge in the quenched case, in the sense that A{p) = 1 
follows by choosing the appropriate nonlocal gauge in the bare vertex approximation and 
that the nonlocal gauge a{k) as a function of the photon momentum k reduces to the 
Landau gauge in the quenched limit A^ ^ 0. Taking into account only the leading and 
next-to-leading order terms in the infrared region, which is the essential region in QED3 
yields a finite critical number of fiavors |TB| 



= 128/37r2 ^ 4.3 



(5) 



which is the same as the leading-order result obtained by Nakatani, and coincides with the 
result obtained by Nash in a different way pSf . In this paper we use the nonlocal gauge as 



derived in [16 



In the present paper the existence of a finite critical number of fiavors is confirmed by 
solving numerically the nonlinear SD equation, using a bare vertex and the nonlocal gauge 
as derived in |T^ in order to keep A{p) = 1, without further approximation. The parity- 
conserving and chiral-symmetry-breaking fermion mass is dynamically generated in QED3 
and there occurs a chiral-symmetry-restoring phase transition at = 128/ (37r^). The order 
of the chiral phase transition is infinite in the sense that the dynamical fermion mass rrid 
and the chiral order parameter (iptp) exhibit an essential singularity at the critical point 



N= Nr 



rrid\ 
a J 



3/2 



exp 



37r 



(6) 



In the presence of CS term, as already reported in |T]], the chiral phase transition turns 
into a first-order transition, in sharp contrast to the ^ = case. This result is obtained by 
both a numerical study of the full SD equation and an analytical study of approximated 
equations using the same scheme as in the absence of CS term. We have also determined 
the critical line A^ = Nc{6) for this transition in the (A^, 6) plane, the two-dimensional 



^The N ^ oo limit (with a being kept finite) corresponds to the weak-coupling limit [3 := 
l/(ee^(e)) — > oo in the lattice gauge theory defined on a lattice with lattice spacing e. 
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phase diagram. In this paper we give the detailed exposition of this result M based on the 



framework of |]16[. The first-order transition implies that the chiral order parameter as well 
as the dynamical fermion mass show a discontinuous change at the critical point on the 
whole critical line extending from {N,6) = {Nc,0) in the phase diagram {N,6). Moreover, 
the critical number of flavors Nc{6) decreases as the CS coefficient 6 increases. 

This paper is organized as follows. In section |I| we give a description of QED3 with a 
CS term. We start from the formulation of QED3 using four- component fermions in the 
reducible representation of the Clifford algebra for the 7-matrices in (2+l)-dimensions. We 
define the chiral and the parity transformations for this theory, and our decomposition of 
the fermion propagator into scalar functions. We also discuss the structure of the gauge 
boson propagator with the nonlocal gauge and give an explicit expression for the vacuum 
polarization in the leading order of expansion. After these preliminaries, we first dis- 
cuss QED without the CS term and subsequently study the effect of the explicit CS term. 
In section we write down the SD equation for the fermion propagator explicitly. The 
SD equation is solved both analytically and numerically. Next, in section we study the 
effect of the explicit CS term by solving the nonlinear SD equation numerically. In section [V| 
detailed comparison of numerical and analytical result is given and the numerical results in 
section ^ are confirmed by the analytical treatment. In order to study the stability of the 
chiral-symmetry-breaking solution, we evaluate the effective potential of Cornwall-Jackiw- 



Tomboulis ||3g] in section 0. The final section is devoted to conclusion and discussion. In 
Appendix A, we give integration formulae which are necessary to perform the angular inte- 
gration to obtain the integration kernel in the SD equation. In Appendix B, we supplement 
details on the calculation of the ultraviolet boundary condition. In Appendix C, it is shown 
that at the stationary point the effective potential obtained in section ^ actually gives the 
SD equation for the fermion mass function in the nonlocal gauge. 

II. QED3 WITH CHERN SIMONS TERM 

In Euclidean space the Lagrangian is 

£ = ^(i ^ + e 4 + me + Tmo)ip + \F^^ - \iet^upA^dyAp + ^gaugo fixing • (7) 

We use four-component spinors for the fermions, and accordingly a four- dimensional repre- 
sentation for the 7-matrices of the Clifford algebra {7^,71^} = —'2,5^^. 



" I ^a^) ' " I -^a') ' " I -^a^) ' 
with o"a(a = 1,2, 3) being the Pauli matrices; furthermore we have the matrix r, defined by 

With such a representation we can define chirality just as in four-dimensional QED, but 
now there are two independent chiral transformations possible, which are defined by the 
4x4 matrices 73 and 75: 
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73 = (^J o) ' ^ 7o7i7273 = (^J • (10) 

Without an explicit mass rrie for the fermions, the Lagrangian is chirally symmetric, but 
the mass term meipip breaks chiral symmetry. Note that the other mass term, rUoipTilj, is 
chirally invariant. 

In this representation, the parity transformation is defined by ip{xo,Xi,X2) — >■ 
Pip{xo, —xi, X2), Af^{xo, xi, X2) {—lY''''-A^{xo, —xi, X2), with P = —^7571. Then the mass 
term moipTip is odd under a parity transformation. Also the CS term is odd under this 
parity transformation, so in the presence of a (bare) CS term parity is always broken, even 
if rrio = 0. The other terms in the Lagrangian, including the chiral-symmetry-breaking mass 
term, are invariant under a parity transformation. 



A. The Fermion Propagator 

The inverse full fermion propagator can be written as 

S-\p) = A,{p) ^ + A,{p)T ^-B,{p) - B,{p)r. (11) 

The functions A{p) and B{p) are scalar functions of the absolute values of the momenta, 
and their bare values are A^ = I, Aq = 0, Bf. = rrie, and Bo = rrio. In order to study the 
fermion propagator nonperturbatively it is useful to define the matrix (projection operator) 

X± = ^(l±r), (12) 
which allows us to rewrite the full propagator as 

S{p) = S+{p)x+ + S-{p)x- 

A+{p)^+B+{p) A^jp) ^+B_{p) 



A\ ip)p^ + Bl (p) Al_ {p)p^ + Bl (p) 
where we have defined 

A± = A,±Ao, (14) 

and 

B± = Be±Bo. (15) 

In this paper we study the dynamical breaking of chiral symmetry, using the SD equation 
for the fermion propagator. Perturbatively this will not happen, but just as in pure QED 
(without CS term), chiral symmetry can be broken dynamically due to nonperturbative 
effects, starting with a chirally symmetric Lagrangian. So in the remainder we have put both 
explicit masses mg and nio equal to zero Q and study the behavior of B± nonperturbatively. 



^We will reconsider this procedure in section VI. 
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Note that in terms of and chiral symmetry means = —B^{p), which gives 

Beip) = 0. 

The order parameter connected with the chiral phase transition is the chiral condensate. 
In the presence of a CS term there are two types of condensates: a parity even condensate 
(ipip) and a parity odd condensate {iprip). Using the decomposition of the propagator in 
terms of x±i 'we can write these condensates as 

(V^V') = (^V^)+ + (^V^)-, (16) 

(^r^) = (^^)+-(V;^)_, (17) 

where 

The general SD equation is given by 

6?k 1 



/(X k 1 
^-Tr[7^^±(fc)r,(p, k)D^,{j> - A;)x±] , (19) 

A±(p2) = l + -y^-Tr[^7;.^±(fc)r,(p,fc)D^,(p-A;)x±]. (20) 



In QED3, the usual truncation scheme for the fermion SD equation is based on the 1/iV 
expansion p6|. The coupling constant has the dimension of mass, and we use the large 
iV limit in such a way that J. and the product N remains fixed: 

e^ = Sa/N, (21) 

with a fixed. In this approximation scheme, the full vertex is replaced by the bare vertex, 
because that is the leading-order contribution in In order to be consistent with the WT 
identity, or at least with the requirement that the vertex renormalization and the fermion 
wave-function renormalization are equal, we use a suitable nonlocal gauge function. 



B. The Gauge Boson Propagator 

We use a Lagrangian with a so called nonlocal gauge-fixing term for the gauge field, 
namely 

C{A- fix - y)) = \Fl - \iee,,pA^d,Ap + i jd^yf{x - y)d,A^{x)d,A,{y) , (22) 

which has certain advantages above the normal gauge- fixing term: it allows for a momentum- 
dependent gauge parameter in the gauge boson propagator. 

It is easy to show that this Lagrangian leads to the following inverse bare photon prop- 
agator 

Dl.-\<D = <l\5,u - ^) + ee,,,qp + a-\q)q,q, , (23) 
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where we have defined the momentum-dependent gauge parameter in the following way 



a-\q) = Jd^xf{x)e-'''^. (24) 

So if we can simplify the actual calculations by choosing a specific momentum-dependent 
gauge a(g), we can justify this afterwards by specifying the gauge-fixing term as 

C = ^ Jd'y fix - y)d,A^{x)d,A,{y) , (25) 

with 

The conventional (local) gauge can be recovered by the choice f{x) = 6{x)/a, as can easily 
be seen from the above formulae. 

We can show that identities relating the different Green's functions, like the WT identity, 
hold with this nonlocal gauge in exactly the same way as with a constant gauge parameter. 
For the (full) photon propagator this implies in momentum space 

qYD^u{q) = -a{q)qu. (27) 

This means that also with this nonlocal gauge the gauge- dependent part of the photon 
propagator is not affected by the interactions, and that the longitudinal part of the vacuum 
polarization is zero. Therefore we can write the full photon propagator as 

D,A<1) = D^{q'){6,^ - ^) + DO{q')e,^,f^ + a{qf-^ , (28) 

q^ \q\ q 

' (g2-n^(g))2+(nO(g)_^|g|)2' ^ ^ 

DO(q^) = (30) 
where 11^ and 11*^ correspond to the decomposition of the vacuum polarization tensor |Q: 

n,.(g) = - ^) n^(g) + 6,.p^n«(g). (3i) 

In the expansion the one-loop vacuum polarization has to be taken into account, 
because this vacuum polarization is of order one: there are fermion loops contributing 
to the vacuum polarization and each loop is of the order ~ Starting with massless 

fermions, e.g. both the parity odd and parity even mass equal to zero, there is no parity 
odd part of the vacuum polarization. The transverse part of the vacuum polarization is 

Uriq) = -a\q\ . (32) 

Therefore the inclusion of this vacuum polarization leads to the following gauge boson prop- 
agator 
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q\\\q\ + a)^ + fc*^^ 



^V) = , (34) 













2 + 02) 



This is the photon propagator we will use in this paper, with a suitable choice for a(g) in 
order to keep the wave-function renormalization equal to one. 



III. DYNAMICAL SYMMETRY BREAKING IN PURE QED 

Without the CS term in the Lagrangian, we have no explicit parity breaking terms, and 
the full fermion propagator will also be parity even: there will be no spontaneous breaking 
of parity ||TO|-|l^ . That means that we only have to deal with one set of two coupled integral 
equations for A{p) and -B(p) 

a f d^k A{k) 



A^p2 J (27r)3 A;2A2(A;) + 52(fc) 
((B.(,)_f^)Mi^ + ^M|^), (35) 



with q = k — p. 



The condition that the wave-function renormalization is equal to one leads to 16 



a{q) = 2q'D^iq) + — - 2 - ^ In ^ . 37 

\q\ g2 \ a J 

The SD equation for the fermions thus reduces to only one nonlinear integral equation for 
the dynamical mass function, which we call m{p). After the angular integration, which can 
be done analytically (see appendix A), the radial integration kernel becomes 

K(p, A;) = + — In f ^±]h±A] 

' max(A; ,p) |A;2 — p2| j^p ^Q,_|_|/j_p|y 

^Infl + ^U , lnfl + ^V (38) 

kp{k-p) \ Oi ) kp{k + p) \ a J 

and the equation for the mass function is 

4a /■°°,, k'^m(k) , , , , , 

^ P = ^ / t^—^I-Mp, k) . 39 
TT^A JO k^ + m'^[k) 

This nonlinear integral equation can be solved numerically without further approximations, 
or it can be solved approximately by using a series expansion for the logarithms. This last 
method makes it possible to convert the integral equation into a second-order differential 
equation, and to study dynamical chiral-symmetry breaking analytically. 
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A. Differential Equation 



In order to solve the equation analytically we make the replacement k"^ + m?{k) — >■ 
k'^ + m^(0), which is a good approximation both for small momenta (where m{k) is almost 
constant), and for relatively large momenta (where both m{kY and m(0)^ are negligible w.r.t. 
fc^). One can also show that this replacement is in agreement with bifurcation analysis, see 
e.g. ||37,25]. For very large momenta {k > a) we neglect with respect to k"^ in the 



denominator completely. 

Another approximation, which is commonly made uses the fact that the integral 
is heavily damped for momenta larger than the mass scale a. All the essential physics 
comes from the infrared part of the nonlinear interal equation. Therefore one uses a series 
expansion for the logarithms in Eq. ( pSf ) for momenta p,k -C a, and introduce a cutoff at 
k = a. However, we are interested in the behavior of the mass function for momenta p > a 
as well, especially when we include the CS term. Therefore we will adopt a slightly different 
approximation, which takes into account the ultraviolet tail of the integral as well. 

For this purpose we expand the integration kernel K(p, k) in powers of 
min(p, k)/ max(p, k). To leading order this gives 

max(p, A;)^ + + Q;max(p, fc) , /i , / tn/ \ 

K{p, k) = 2 \.3. ^ , l\. ^-TTi In (1 + max(p, k)/a) . (40) 

max(p, /cj'^(a + max(p, fcjj max(p, /cj^ ^ ^ 

Differentiating Eq. (^) with this kernel leads to a second-order differential equation 

f{p)m'\p) + g{p)m\p) = , (41) 

with 

and two boundary conditions, infrared and ultraviolet ones, just as in the usual approxima- 
tion scheme. It is easy to show that the behavior of this differential equation in the infrared 
region is exactly the same as one would obtain by expanding / and g for small momenta 
directly. 

For momenta j9 3> a we expand the functions f{p) and g{p), which gives to leading order 
in p 

IGof 

p^m"{p) + 3p'^m'{p) + -^—m{p) = . (44) 
To solve this we substitute a solution 

oo 
1=0 

and again taking into account leading order only we find a = or a = —2. It is easy to 
see that a = is the solution corresponding to explicit chiral-symmetry breaking, since 
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this leads to m{p) — > mo 7^ in the ultraviolet region. Without a bare mass, the correct 
solution is a = —2, which shows that the dynamical mass function falls off very rapidly in 
the ultraviolet region. This is also consistent with the ultraviolet boundary condition, but 
it does not provide a normalization condition. 

With this knowledge we now consider integral equation in the infrared region, without 
neglecting the ultraviolet tail of the integral: we use 

m(p) = a^^TT^ (46) 

for p > a, based on the ultraviolet behavior we have just found, and normalized in such a 
way that the solution is continuous aX p = a (we do not require the continuity of derivative, 
m'{p), ai p = a). Using the usual approximation 0] for p < a, we arrive at the integral 
equation 

, , 32 Pm{k) 1 8a m{k) 

Sn'^Jy Jo + m^(0) max(A;,p) ti'^M Ja 

which finally reduces to 

32 k'^m{k) 1 8m{a) . 

""^^^ " 3^ io F + m2(0) max(fc,p) ^ ' ^ ^ 

This leads to exactly the same second-order differential equation in the infrared region 

p'm^p) + 2pm' {p) + / = , 49 

as that obtained by neglecting the ultraviolet tail completely, with also the same infrared 
boundary condition 

m'(0) = . (50) 
The general solution of Eq. (^9]), satisfying the infrared boundary condition, is given by 

m{p) = m(0) 2Fi(a+, a_, | -pV^(O)') , (51) 



with a-i- = i ± ^iy Nc/N — 1 and a critical number of fermion flavors 

iVc = ^ , (52) 

above which there is no chiral-symmetry breaking. The only difference is the "ultraviolet" 
boundary condition at p = a. This boundary condition now becomes 

, , ,, , Smia) 

m{a) + am (a) = ^^^^ , (53) 

due to the inclusion of the ultraviolet tail of the integral. This boundary condition leads to a 
slightly different normalization of the infrared mass m(0) (see appendix B for more details) 
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m(0) 



a 



exp 



-2n + 2(f) 
NJN - 1 



(54) 



with 



arg 



T{l + -^NjN-l)T{a_y (^l--a. 



(55) 



Close to the critical number of flavors this can be expanded to give 

-271 



m{0)/a = exp 



.JNJN- 



1 2 
= + 31n2 + -7r + - 
1 2 7 



(56) 



which is almost the same as if one neglects the ultraviolet tail, in which case the term | in 
the exponent is absent. 

However, this is not the only solution for the boundary condition; it is known that there 
are infinitely many solutions if < A^^c which behave in general as 



m{0)/a = exp 



-2ri7r , 1 2 

+ 31n2 + -7r + - 

NJN - 1 ^ 



(57) 



close to the critical number of fermion flavors. The solution with the largest value ofm/a < 1 
corresponds to the ground state, since this has the lowest energy. That means n = 1 in the 
above equation, and one can also show that the other solutions have an oscillating behavior 
at large momenta; only the solution with n = 1 is a nodeless solution. 



B. Integral Equation 

Alternatively, we can solve the nonlinear integral equation for the mass function, Eq. (p9[), 
numerically, without further approximations. The integral can be calculated numerically 
without any cutoff, since the integrand falls off rapidly for large momenta and the integral 
is finite. Solving the integral equation iteratively leads to qualitatively the same result as 
the previous analysis. 

There is a critical number of fermion flavors below which there is a chiral-symmetry- 
breaking solution. The critical number is the same as we have just found analytically 



Nr 



128 
3^ 



4.32, 



(58) 



and also the behavior of the infrared mass m(0) is similar as before. In Fig. |l] we have 



plotted the infrared mass, or actually \JNJN — 1 ln(m(0)/a), versus \JNJN — 1 for the 
full integral equation together with the analytical result. We can see that in both cases the 
mass behaves like 



m(0)/a = exp 



-2cTi 



NJN 



+ (T2 



(59) 
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the only difference is the value of the constants (jj. The analytical result is (Xi = vr and 
(J 2 = 3.94, whereas the full nonlinear equation gives di = 3.1 and (T2 = 1.8. So the main 
difference between the analytical solution of the approximated equation and the numerical 
solution of the full equation is an overall normalization factor. This difference is due to the 
fact that in the first case we have made several approximations, but these approximations 
turn out not to be essential for the behavior of the infrared mass near the critical point. 
Both in the infrared and in the ultraviolet region the numerical and analytical solutions have 
the same behavior, see Fig. ^ Of course, these approximations do make a difference in the 
region where p is of order a, but this is to be expected since we expand to leading order 
only in p/a and a/p respectively. 



C. Chiral Condensate 



The order parameter of the chiral phase transition is the chiral condensate 



TT^ JO 



Cl/b ~ 



^2 _^ rn'^(^k) 



(60) 
), which 



In order to get an explicit analytical formula for the condensate, we can use Eq. 
allows us to write the condensate in terms of the mass function aX p = a. Using the same 
linearization as before we get 



4 



32 



amial 



(61) 



From this equation we can calculate how the chiral condensate behaves close to the critical 
number of fermion flavors, see appendix B, 



{ipip) _ 608 r(l/2) 
~ 77r2r(l/4)2 



exp 



-3n 



NJN 



9, o 3 3 
= + -ln2 + -7r + - 
1 2 4 7 



(62) 



Alternatively, we can calculate the condensate numerically, using the numerical solution of 
the full nonlinear integral equation, which shows a similar behavior, see Fig. |I[ 



IV. EXPLICIT CHERN SIMONS TERM 



In this section we add a CS term for the gauge field to the Lagrangian. This breaks 
the parity explicitly, and gives rise to a parity odd mass term for the fermions, as well as 
a parity odd part of the gauge boson propagator. We use the gauge boson propagator as 
discussed in section [11 B| , with the inclusion of the leading-order vacuum polarization. 

Again we will use a nonlocal gauge-fixing term. For the momentum-dependent gauge 



function we use flG 



a{q) = 2q^D'^{q) 




1 H — arctan 



e\q\ 



+ a\q\ + 6"^ 



(63) 
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which leads to a wave-function rcnormahzation almost equal to one, the deviation of one is 
proportional to 9B±, which can be kept very small. 



A. Schwinger— Dyson Equation 

With the above gauge boson propagator and gauge function, the SD equation can be 
rewritten into two decoupled sets of two coupled integral equations for and S+, respec- 
tively for A_ and This leads to the following set of equations: 

8a r d^k 1 

J {27r)n^ Al{k) + BUk) 

(^B^{k){2D^{q) + 2A^ik)D^{q)^'j . (65) 

The even and odd parts of the scalar function B{p) are 

B, {p) = {B+{p)+B_{p))/2, (66) 

and 

Bo{p) = {B+{p)-B4p))/2, (67) 

in terms of B±{p), and similar for A{p). In analyzing these equations it is important to 
observe that once we have found a solution for A^{p) and B+{p), we automatically have also 
a solution for A_(p) and B_{p): namely the set A-{p) — A^{p) and B_{p) = —B^{p). That 
means that we can always construct a chirally symmetric (but parity odd) solution of the 
SD equation, with B^{p) = 0. The question of dynamical chiral-symmetry breaking turns 
into the question whether or not there exist two (or more) solutions of the set of integral 
equations. 

Without the CS term there is dynamical chiral-symmetry breaking only for N < Nc — 
128/(37r^) as we have just seen. We therefore expect a similar situation in the presence of 
the CS term, at least if the parameter 9 is small. That means that for > N^, we have 
only the chirally symmetric solution of the above equations, but for N < we expect that 
there are (at least) two solutions for both B^ and i?_ possible, which can be distinguished 
by their behavior under the chiral and parity transformations and by their behavior in the 
limit ^ i 0. 

It is also essential to note that in the presence of an explicit CS term for the gauge field 
in the Lagrangian there is no trivial solution B = 0, as there would be without the explicit 
CS term. Due to the explicit breaking of parity, the fermions always acquire a parity-odd 
mass Bo, even if the explicit odd mass rrio is zero. At large momenta the CS term in the 
SD equation dominates (at least if there is no explicit mass term present), which leads to 
an ultraviolet behavior of the mass functions 

B+{p)r^-B_{p)r^^-, (68) 
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at a perturbative level, whereas in the absence of the CS term the mass function m{p) falls 
off like as we have seen in the previous section. We will discuss this point in more 

detail when we are studying the breakdown of chiral symmetry analytically. 



B. Numerical Results 

We have first solved the set of coupled integral equations Eq. ( p^ ) and ( |B5D numerically 
without further approximations. Depending on the values of the number of fermion flavors 
and the CS coefficient 6, there is only one (chirally symmetric) solution in most parameter 
space, but we found two solutions, allowing for a chiral-symmetry-breaking solution, for 
small 6 and small A^, in agreement with the above expectations. Using the notation B±{p) 
and B±{p) for the two different solutions of the above equations, we have found the following 
solutions: 

1. 5+(0) = -5_(0) = C(m(0)) > 0, with B+{p) = -B_{p) > m{p), 
B^{p) = and ^^(O) = C(m(0)) 

for < N(.{9 = 0) we find that B±{p) and thus Bo{p) go towards the nontrivial 
solution m{p) in the limit 10; 

for > Nc{6 = 0) we find that B^{p) = Bo{p) = 0{6), which vanishes in the limit 

2. 5+(0) ~ 5_(0) = O(m(0)), with B+{p) > m{p) > B_{p) 
B^{p) = 0{m{p)) and Bo{p) = 0{9) 

This solution exists for values of A^ < Nc{6 = 0) and small (compared to a) values 
of 9 only. There is a critical number Nc{6) < Nc{0) for given 6, or critical 6c{N) for 
given A^ < Nc{0). In the limit 6' | both B± and thus Bg go towards the nontrivial 
solution m{p): B^{p) | m{p) and B_{p) f m{p) in this limit. The parity breaking 
solution Bo{p) vanishes in this limit. 



3. 5+(0) = -5_(0) = -O(m(0)) < 
The chirally symmetric combination, 
values of 6 only. 

4. S+(0) ~ 5_(0) = -O(m(0)) 

The fourth possible combination of -B+ 
small values of 9 only. 



which exists for small values of A^ and small 



and -B_ , also existing for small values of A^ and 



The typical behavior of the numerical solutions for A and B is shown in Fig. ^ As 
expected, A± is indeed very close to one, due to our choice of the gauge function. The 
solution Bj-{p) exists for all values of both A^ and 6, allowing only a chirally symmetric 
solution. The other solution, B±{p), is the interesting one, leading to dynamical chiral- 
symmetry breaking. The iterative process of solving the integral equation numerically does 
not converge to a (second) stable solution B± for all values of A^ and 6: we can find this 
solution only for A^ below some critical value (depending on 6) and 6 below some critical 
value (depending on A^), thus showing a chiral phase transition at some critical 6 and A^. 

Our numerical results all indicate strongly that this phase transition is a first- order phase 
transition, in contrast to the infinite-order phase transition in pure QED (^ = 0), although it 
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is numerically very difficult to establish the type of phase transition at the critical values of 
N and 6. For fixed 6 and small A^, there exists a second solution, B±{p), as can be seen from 
Fig. ^. The value at the origin -B-(O) decreases rapidly for increasing and this solution 
"disappears" (without B±{0) going to zero) at some critical value N < Nc{9 = 0) , and also 
Be{0) does not go to zero at this critical value. Considering Be{0) as the order parameter 
for the chiral phase transition, this corresponds to a first-order phase transition, in contrast 
to the infinite-order phase transition at ^ = 0. 

If we look at the behavior of B_{0) and -Be(O) at fixed N < Nc{6 = 0) and increase 6, we 
see a similar situation. For (very) small values of 6 we find the two solutions B and B leading 
to chiral-symmetry breaking. For increasing 6 we find that -B(O) decreases, and disappears 
at some critical value 6c{N), without going to zero at this critical value, see Fig. ^ 



C. Condensate 

Although both functions B±{p) behave like 6/p in the ultraviolet region, leading to 
divergent integrals for the condensates {ipi/^)±, the chiral condensate is convergent due to 
the fact that the leading-order contributions in B+{p) and B_{p) (or B_{p)) cancel. So the 
chiral order parameter is 

,j , _ 1 ( PB+jk) k^B4k) 

~ 7r2 Jo \Al{k)k^ + Bl{k) + i2 (^)p + Bl{k) 

The chirally symmetric combination, {B^{p), B^{p)) with B^ = 0, gives (ipip) = 0, as would 
be expected from a chirally symmetric solution. Note that the other condensate, (iprip), 
is actually logarithmic divergent because the leading-order contributions add up. Once we 
have the numerical solutions, it is straightforward to calculate this chiral condensate as well, 
see Fig. ^. The behavior of the chiral condensate also indicates a first-order phase transition. 

However, as we mentioned before, it is numerically very difficult to establish such a 
ffist-order phase transition. In order to conffim that it is indeed a ffist-order transition 
and to determine the critical parameters Nc and 9c, we have to study the phase transition 
analytically. We can do this by solving the integral equations analytically, after some further 
approximations analogously to the approximations leading to analytic solution in the pure 
QED case. 




V. ANALYTICAL STUDY 

In order to see whether there is indeed a ffist-order phase transition, we have analyzed the 
SD equation analytically, after some more approximations. Based on the fact that A{p) = 1 
exactly if 6 is zero, and very close to one for small values of 6 (small compared to a), we 
replace A{p) by one, so we get the following equation for B± 

d^k 1 



B±{P) = ^ j (^2Tifk^^Bl{k) 



B^{q)[2D^{q) + ^)T2DOiq)^\ . (70) 
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The error we make in neglecting the effects of the wave-function renormahzation (even 
if we use the nonlocal gauge) is of the order of 9B±, which is (close to the critical number 
of fermion flavors in the absence of the CS term) of the order of 9"^, because the parity odd 
mass, generated by the CS term, is of order 9. So for a consistent approximation in order to 
get an analytical solution of the equation, it is enough to expand the kernel in 9 and neglect 
all terms of order 9'^ and higher. Expanding the nonlocal gauge in 9 gives 

a(g) = ^ + ^- 2- i^ln (^) + 0{9^) . (71) 

So to order 9 it is the same as the nonlocal gauge without the CS term, see Eq. (0). By 
inspection of the transverse part D'^[q^) of the photon propagator it is easy to see that up 
to order 9 this is also exactly the same as without the CS term. That means that the second 
term on the RHS of Eq. ([70|), proportional to D^{q), is the only 0{9) contribution in the SD 
equation for B±, and that for the other terms we can use the same kind of approximations 
as in section ^ on pure QED. 

So in order to study the problem of dynamical chiral-symmetry breaking analytically, we 
replace k"^ + B^{p) by k"^ + M^, where M± = B±{0), in the denominator of the integrand. 
In general we get 



B^(p) = r dk 'f.^^Kip, k) 

^^^^ Jo k^ + Ml ' 



Aa9 f°° k'^dk k ■ {k - p) 



Ti^N Jo k^ + MlJ-i {k-py{\k-p\+aY' 

with the kernel given by Eq. (|38D , and z = cos(f), with the angle between the vectors p 
and k: p ■ k = pk cos0. For momenta p < a this reduces to 

32 k^ B±{k) 8a B±{k) 



Stt^A^Vo P + M| max(p. A;) vr^A^ 7a P 

Aa9 /•°° k'^dk /"i k ■ (k - p) 

dz— — ^— ^. (73) 



Jo k^ + MlJ-i {k-pf{\k-p\+ay 



In order to get this equation, we have used the same approximations as in section ^ for the 
first term of the integral, which is independent of 9. 

A. Explicit Chern Simons Contribution 

The last term on the RHS in Eq. ([7^), proportional to 9, 

4a6' /-"^ k'^dk k ■ {k - p) , , 

^HP) = ^y^ ^2 + Ml J -I {k - py{\k -p\ + a)2 ' ^ ' 



^This should be compared with the analysis by Hong and Park [40| where the order 9 term was 
neglected from the beginning. 
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can be calculated analytically, see appendix A. For the radial integration we expand the 
integration kernel for p < q and p > q and take into account only the leading-order terms in 
mm{p, q)/ max(p, q). Furthermore, we neglect M± with respect to a, which is justified close 
to the critical number of fermion fiavors and for small values of 6. So in the far infrared 
region the leading contribution coming from the explicit CS term behaves like 

le{p) = :^e + 0{p), (75) 
and in the far ultraviolet region the leading-order behavior is 

This means that in the ultraviolet region the CS term will dominate, since we know that the 
mass function B{p) in the absence of the CS term behaves like in the far ultraviolet. 
Higher order contributions in mm{p, q)/ max(p, q) will change this result only quantitatively, 
but not affect the general behavior. 

Based on these expansions, we use an analytical (continuous) approximation for the 
contribution from the 9 term 

for T) Oi 
3^f^-^l for«>a 

for small values of 9 compared to a. 



B. Ultraviolet Behavior 

The next thing we have to calculate before we solve the integral equation in the infrared 
regoin, is the ultraviolet tail 

^~dt%P. (78) 



Tl^N J a 

but for this purpose we have to know the ultraviolet behavior of the mass function. From 
the integral equation we can derive a second-order differential equation for momenta p ^ a 

miB'lip) T im + 9{P){B'^{P) T I'eip)) = , (79) 



with f{p) and g{p) given by Eq. (^21) and ( ^3]) respectively. In the ultraviolet region, these 
functions behave like 

3j^^ 

fip) = — T' 9ip) 



4 ' -^-^ 4 ' 

I'eip)-^, I'M 



The leading ultraviolet behavior of the solution of this equation is either constant (which 
would correspond to an explicit mass term in the original integral equation), or 
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B,(p)^±IM^^%, (80) 

which is the correct solution in this case, and also in agreement with the ultraviolet boundary 
condition. This is the same as the perturbatively dominant behavior, as could also be seen 
directly from the original integral equation, assuming that the mass function falls off rapidly 
enough in the ultraviolet region. In order to have a continuous solution at p = a, we use 
the next-to-leading order term, so we have 

^ / X 32 a ^ / / X 320 \ 

for p > a. 



C. Analytical Results 

So we arrive at the integral equation for momenta p < a 

„ , , 32 A;2 BJk) 8 f 169 , A 80 

= 3^ io P + M^max(p,fc) + (^MV + ^^^^V ± ^ ' ^^'^ 

which we can now solve analytically by converting it to a second-order differential equation 
with boundary conditions. 

The second-order differential equation is the same as without the CS term, Eq. (^), and 
also the infrared boundary condition is the same, the only (but essential!) difference is the 
"ultraviolet" boundary condition at p = a 

+ + (83) 



which should be compared with Eq. (53). It is also important to keep in mind that the 
normalization condition is 

5±(0)=M±, (84) 

and that the SD equation does determine the sign of the mass function as well. This is 
obviously not the case if = 0. 

The general solution of the differential equation satisfying the infrared boundary condi- 
tion and the normalization condition is 

B±{p) = M± 2i^i(a+, a_, | -p^Ml) , (85) 



where a± = j{l ± i^jNc/N — 1). The ultraviolet boundary condition leads to the condition 

18 3 

M± (a+, a_, -; -«Vm|) - 2^1(0+, a_, -; -«Vm|) 
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In order to determine M, we can plot the LHS of the above equation divided by 
(^1 + g^T^^ as function of M for a given value of A^, which gives us automatically 9 
as function of M, see Fig. |^a. From this figure we can see that there are three solutions 
possible for and B_. Two of them correspond to the two solutions which we have found 
numerically, and the third is an oscillating solution which in the limit of vanishing 6 cor- 
responds to the first oscillating solution in pure QED, = 2 in Eq. (^). For extremely 
small values of M and 6 there are more oscillating solutions, in the region around the origin, 
but these are numerically unstable, and correspond to higher excites states, just as in pure 
QED. 

Using the same notation as for the numerical solution, we have -B+(0) > and -B+(0) < 0, 
and for B^ the opposite signs. The set is the chirally symmetric solution which is 

present for all values of ( A^, 6) . The chiral-symmetry-breaking solution is given by the set 
(-B+, B_) (or vice versa, the set (-B+, B_)). In the limit 6^0, this chiral-symmetry-breaking 
solution behaves in the following way 

Be^m, (87) 
Bo = 0{e) . (88) 

If we increase 6 from zero at fixed A^, then the absolute value of the second solution, B±, 
decreases (whereas the other solution increases), until at a critical value 6c this solution 
coincides with a third solution, and disappears, but does not become zero at the critical 
point, see Fig. |a. This clearly signals a first-order phase transition, as was also suggested 
by the numerical results. For comparison, we also showed our numerical data in the same 
figure, where we use a different value of a for our numerical and analytical calculation, anum 
and ttan respectively, in order to have equal numerical and analytical values of m{0) / annm 
in the absence of the CS term (note that the axes are 9/annm and |-B(0)|/Q;num respectively). 
This shows clearly that both our numerical and our analytical results are in good agreement 
with each other. 

We can also plot M versus for some fixed values of 6, see Fig. ^3. From that figure we 
can see that if we increase A^ for fixed 6, the chiral-symmetry-breaking solutions disappear 
if A^ exceeds some critical value Nc, which decreases rapidly as a function of 6. This figure 
shows that the chiral phase transition is first order in this direction as well: increasing A^ 
beyond Nc{6) makes the second (and third) solution disappear, but at the phase transition 
neither B nor B^. = [B^ + B_)/2 (which can be regarded as the order parameter of the 
chiral phase transition) becomes zero. In this figure we can also see that in the limit ^ 
the critical value goes towards A'c(O) ~ 4.32. 

The critical parameters and 9^ can be calculated by differentiating the ultraviolet 
boundary condition, Eq. (p6D, with respect to M. This leads to an equation for as 
function of N^. (remember that M does not vanish at A'c) 

2Fi(a+, a_, -a" /Ml) + 2^1(0+, a_, ^; -a" /Ml) 

16 3 

2Fi(a+, a_, -; -a" /Ml) = , (89) 



37im " ^' 2' 

which can be used as input for the ultraviolet boundary condition itself in order to calculate 
6c- Although we do not have an explicit form for 6c as function of A^, we can calculate it 
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numerically, and have shown the critical line in Fig. |^. In this figure we also show some 
estimates of the critical parameters based on our numerical calculation. This shows that the 
numerical and analytical results are qualitatively in good agreement with each other, and 
the only difference is an overall scale factor (just as in pure QED). 

For small values of Mc and Nc{9) close to Nc{0) ^ 4.32, we can expand the above equation 



in JNc{0)/N - 1, leading to 



Mc/ a = exp 



-2 IT 



1 



7r + 3 ln2 



8 
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(90) 



/iV,(0)/iV-l 2 

(see appendix B), which can be inserted into the series expansion for 6. To leading order in 
A''c(0)/A^ — 1, this gives the following expression for 9c 



4481(1/2) 

He — : — r-T77 exp 

' 75r(l/4)2 ^ 



-3tt 



Nc/N 



3 9, 4 

= + -77 + - In 2 

1 4 2 7 



(91) 



As mentioned before, we can also calculate the chiral condensate which is well-defined, 
even without cutoff. The chiral condensate is 



_ 1 roo 



77^ Jo 



dA; 



e B+{k) 



+ — 



Alik)k^ + BUk) ' Al{k)k^ + B^_ik) , 
In the approximations we are using here, this reduces to 



77^ Jo 



" f k^B+jk) e B_{ky 

[k^ + Ml ^ k^ + Ml 



dk{B+{k) + B^{k)) 



(92) 



(93) 



which can be calculated analytically, using Eq. (|8^) at p = a and the ultraviolet expansion 
Eq. M for 5±. This leads to 



i^ij) = -^(B+{a) + B4a)) { + - 



77' 



32 



4 



(94) 



with -B+(a) and B_{a) determined through Eqs. (q^) and (|86D . In Fig. |] we show the 
chiral condensate (together with the -Be(O)) as obtained both numerically and analytically 
for 6* = and 6 = 10~^ as a function of A^. Again, they are qualitatively in good agreement 
with each other, and show a discontinuity in both the chiral condensate and Bc{0) at the 
critical point. 

In Fig. ^ we show both Bc{p) and Bo{p), as obtained numerically and analytically. The 
numerical and the analytical solutions have the same behavior in both the infrared and the 
ultraviolet region; around p = a there is of course a kink in the analytical solution due to 
the approximations we have made. There is a scale difference between the analytical and the 
numerical solution for the chiral-symmetry-breaking solution, just as in pure QED; in fact, 
if we compare this figure with Fig. ^ we can see that the parity-even mass function is almost 
the same as the dynamical mass function in pure QED. Furthermore, this figure shows very 
well the difference in the ultraviolet between the even and the odd mass function: the odd 
mass function behaves like 6/p, whereas the even mass function Bc{p) = {B^{p) + B^{p))/2 
behaves like 1/p^, due to the cancellation of the terms proportional to ±6 /p. 
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VI. CJT EFFECTIVE POTENTIAL 



We consider the effective action of Cornwall-Jackiw-Tomboulis PB[ (CJT) which is given 



by 



with 



Tit 



r[5] = ro[5] + ri 



-Tr[Ln[5-^5o] + ^0-^5- 1] 

,2 



(95) 

(96) 
(97) 



where is the full photon propagator, and Sq and 5* are the bare and the full fermion 
propagator respectively. At the stationary point 6T[S]/SS = 0, this effective action actually 
gives the SD equation: S^^ — Sq'^ = —e^Ti['y^D^,y'y,yS], as is shown explicitly in appendix 
C. Note that we have chosen the bare vertex to write the effective action, which can be 
justified by taking the nonlocal gauge. Furthermore we remark that the effective action is 
normalized so that To[S = So] = 0. For a review of the effective potential, see e.g. 

We can rewrite the above expression in terms of the fermion propagators S+ and S- 



S = S+X+ + S-X- ■ 



(98) 



By using the fact that x+ ^i-nd X- are projection operators and the relation Tr Ln = Ln Det, 
it is not difficult to show that both To[S] and ri[5'] are decomposed into two parts 



0,1 [ 



ro,i[s+] + ro,i[s_] 



where 



and 



TolSi] = -TT[Ln[S^'So] + S^'S^ - 1] 



(27 



(2vr) 



(99) 



(100) 



;ioi) 



Now we define the effective potential V[S] by dividing the effective action by the space- 
time volume / (Px: V[S] = T[S]/ J dPx. In momentum space, ^[5-1-] is given by 



N 

2^ Jo 



p^dp 



In 



'p'Ai{p) + Bi{py 



+2 



p'A^{p)[A^{p)-l]+Bl{p) 



p^aUp) + bUp) 

For Vif^i], after calculating the trace in the integrand, we arrive at the result 



(102) 
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2D 



(27r)3 J (27r)3 + BUp)WAl{k) + BUk)] 

+ ^(2(, .p)i,.k)-p.k)^ A,ip)AAk) 

-{2DT{q) + ^)5±(p)5±(A;) ± £^[q . pA^{p)B^{k) - g • kB^{p)A^{k)]^ , (103) 

with q = k — p. At this stage we adopt the nonlocal gauge a{q) which leads to no wave- 
function renormalization for the fermion: A±{p) = 1. Then we can write 



V[B] = V[S] 



A±=l 



Vo[B+] + Vo[B_] + V,[B^] + V,[B. 



where 



Vo[B, 



N 

2^ Jo 



p^dp 



In 



>2 + 



+ 2 



Blip) 



p2 + Blip) 



and 



Vi[B, 



d^p f di^k 



X • 



(27r)3 7 {2ny [p^ + Bl{p)][k^ + BUk)] 
-2Dr{q) - #1 B^{p)B^{k) ± ^^[(g ■ p)B^{k) - {q ■ k)B^ 



(104) 



(105) 



(106) 



Actually this reproduces the SD equation for B± in the nonlocal gauge at the stationary 
point 5V[B]/ 5B± = as shown in appendix C. 

At the stationary point of B± satisfies the SD equation and hence B± is given by 

the solution to the SD equation: B± = 5^°'. Therefore, at the stationary point, we obtain 



Vi[B 



d^p B±{p)[B±{p) -m±] 



(27r)3 p^ + BUp) 
Hence the effective potential at the stationary point is given by 



(107) 



J soil 



V[B,,i] = V[Bf] + V[B'_! 



V[B, 



d^p 
(2^ 



In 1 



Bl{p)\ B±{p)[B±{p)+m±] 



p2 + Blip) 



(108) 



Note that the function g{x) = ln(l + x) — is positive and monotonically increasing in 
x{> 0). Therefore, when m± = 0, 



V[B. 



1 



soil 



27r2 Jo 



p^dp 



^9 



(109) 



is non-positive, ^[-8^0/] ^ 0, and all the nontrivial solutions have lower energy than the trivial 
ones (in the case of pure QED). Therefore we can determine the ground state, namely the 
solution which minimizes the effective potential, by using Eq. (|109|) . 
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Without actually calculating the integral, we can easily conclude that in the presence 
of an explicit CS term, the chirally symmetric solution {B^,B_) gives a lower effective 
potential and is thus favored above the chiral-symmetry-breaking solution {B^,B_). Since 
the effective potential is the sum of two non-positive integrals, each of which depends on 
B^ or B^, this effective potential is minimized by a symmetric set of solutions and 
not by {B^,B_). So in the presence of an explicit CS term, the chirally symmetric phase 
is the ground state, even if the number of fermion flavors N and the CS parameter 6 are 
below their critical values for dynamical chiral-symmetry breaking. This is quite surprising, 
since it is well-known that without the explicit CS term the chiral symmetry is broken for 
N below the critical number, as can also be seen from the effective potential, Eq. ( |109|) . 

In solving the SD equation we have put the external source rrif. and rrio to zero from 
the beginning. As a result, the solution for the homogeneous SD equation has no specific 
direction for the dynamical mass to be generated. Actually, if a solution is found for then 
another solution i?_ = — -B+ is automatically obtained, and it is this solution which gives 
the lowest effective potential. In order to study the spontaneous breaking of a symmetry in 
general, one can introduce an external source JO which breaks the symmetry in question and 
subsequently consider the limit of removing the external source. If the symmetry is broken 
even in this limit, which is signaled by the non-vanishing order parameter = \imj^Q{0)j, 
then it is said that the symmetry is spontaneously broken. In taking this limit one must 
specify from which direction the external source is decreased to zero, J 0. 

In this case we must consider the two limits: rrie —>■ 0, rrio — * 0. For simplicity we keep 
rrie and rrio positive, rrie, rrio > and consider the limit: rrie | 0, | 0, without loss of 
generality. Then there are five cases to be examined, which are given below. 

explicit masses condensates solutions of SD eq. 

me>mo>0 m+>m_>0 {ip'^) + > (tptp) - > 5+(0) > 0, S_(0) > 

me>0,mo = m+ = m_>0 + = {ip^) S+(0) > 0, 5_(0) > 

mo > me > m+ > -m_ > > >_{ipip)- S+(0) > 0, B_{0) < 

mo>0, me = m+ = -m_ > (#)+ = -(V^^)- > 5+(0) = -5_(0) > 

mo = me> m+ > 0, m_ = {^^)+ > 0, -B+(0) > 0, 

(ipip)- undetermined both -B-(O) and -B-(O) 

Therefore the two solutions {B^,B^) and {B^,B^) are realized in two different limits: the 
limit m+ = — m_ J, leads to the solution {B+, B_), whereas the limit m+ = m_ J, yields 
the solution (5+, B_). This last set of solutions is of course only possible if the parameters 
N, 9 are in the chirally broken phase. 

From this consideration, the two sets of solutions cannot be realized simultaneously. 
The result which solution is realized depends on the ordering of taking the vanishing limit 
keeping the relation between rrie and rrio. So in the presence of the CS term, taking the 
limit rrie > i 0, we will get the solution {B^,BJ) which has a higher value for the 
CJT effective potential than that of the "trivial" solution That means that this 

solution is a quasi-stable solution of the SD equation, and will eventually decay into the 
"trivial" solution. This quasi-stable solution does no longer exist beyond a certain value 
of 6', the critical value 6o{N) depending on N. At this point the system jumps directly 
to the state expressed by the solution {B^,B^), which is energetically more favored, in a 
discontinuous way. 
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VII. CONCLUSION AND DISCUSSION 



In (2+l)-dimensional QED with flavors of four-component Dirac fermions, we liave 
solved the SD equation for the fermion propagator in the nonlocal gauge both analytically 
and numerically. In the absence of the bare CS term, we have shown the existence of a 
finite critical number of flavors Nc = 4.3, below which the chiral symmetry is spontaneously 



broken, in agreement with previous analyses |l27| , ^ , p!6| . In the presence of CS term, 9 



we have obtained the critical line extending from the critical point (A^, 6) = {Nc, 0) in the 
two-dimensional phase diagram {N,6). Here a quite remarkable point is that, no matter 
how small 6 is, this phase transition turns into a first-order transition when 6' 7^ 0, although 
the critical point {Nc, 0) in the absence of CS term is a continuous (infinite-order) phase 
transition point. Therefore the point {Nc, 0) is the only one point on the critical line which 
exhibits continuous phase transition. 

We have shown that the chiral-symmetric solution {B^,B_) and the chiral-symmetry- 
breaking one B_) of the SD equation are stationary points of the effective potential of 
the CJT type. Here it should be remarked that the trivial solution {B^ = 0, B_ = 0) as 
a trivial stationary point (equivalently the trivial solution of the SD equation) is possible 
only when ^ = 0. Although the symmetric solution {Bj^,BJ) gives lower effective potential 
than the symmetry-breaking solution {B^,B_), we have shown that it is possible to take 
the limit of removing the external source, me, mo | 0, so that the spontaneously chiral- 
symmetry-breaking solution {B^,BJ) is realized. However it is not yet clear whether they 
give local minima, local maxima, or possibly a saddle point of the effective potential and 
whether or not (i?+,i?_) gives the absolute minimum. In order to confirm this issue, it is 
necessary to calculate the second functional derivative: 

Vcjt[B+,B_] (110) 



5B^{p)5B^{q) 



at the respective stationary point The system at an excited state {Bj^,BJ) might be 
quasi-stable and might decay into the more stable state (-B+, BJ) in a finite time interval. 
The stability of the stationary point will be discussed in more detail in a subsequent paper. 

It is interesting to see our result from the viewpoint of the lattice gauge theory where 
the lattice spacing a corresponds to the inverse of the UV cutoff A, a ~ 1/A. It is well 



known |^4[ that the nontrivial continuum limit can be possibly taken only at the second- 
and higher-order transition point where the correlation length diverges (in units of the 
lattice spacing). Hence our result implies that the meaningful continuum limit from the 
broken phase in the nonperturbatively regularized three-dimensional gauge field theory can 
be taken only at the point {Nc, 0), i.e. in the absence of the bare CS term. This is a novel 
feature of the three-dimensional gauge field theory with a bare CS term, which seems to 
be overlooked so far. Our result may have important implications in the application of the 
three-dimensional gauge theory in condensed-matter physics as a long-wavelength effective 
theory of the microscopic model. This will be discussed in a forthcoming paper. 

This result is obtained by using a consistent expansion of the full gauge boson propagator 
and the vertex in 1/A^. In order to satisfy the WT identity (up to terms of order 1/A^ and 
of order of the dynamically generated mass), we have adopted the nonlocal gauge function 
which allows us to neglect effect from the wave-function renormalization. A more accurate 



25 



approximation would include a full calculation, similar to the analysis in in pure 

QED. Since in pure QED these corrections do not change the result qualitatively, but 
only quantitatively (leading to a slightly different critical number of fermion flavors), we do 
not expect that those corrections will change our present result essentially. The question 
of gauge covariance can (in principle) be recovered by applying the Landau-Khalatnikov 
transformation rules to the various Green's functions Therefore we expect our present 
results to hold also in a more elaborate approximation of the SD equation. 
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APPENDIX A: INTEGRATIONS 



1. Angular Integration Formulae 



First we note the following relation 
d0sin0/(g) = 



1 \ rp+k 

dzf{q) = — / qdqf{q) , 

1 pk J\p-k\ 



(Al) 



where = {p — k)^ = p'^ + k'^ — 2pkz with z = cos 0. Then we obtain the following formulae 

2 min {p, k) 



-1 q 

-1 pk \p — k\ 

11 1/1 1 

dz 



pk ' 

1 , p+k 
In ■ 



1 q^ pk \\p — k\ p + k) ^ 



' dz\ In f 1 + ^ 

-1 q^ \ a 



1 

pk 
+ 



2a \ \p — k\ p + k J 

1 ( \p — k\ 

In 1 + ' 



2«2 



In- 



1 + 



\p~k\ 



1 + 



p+k 



dz- 



1 



2|p- 

1 , p+k + a 
In ■ 



a 



2{p + ky 



a 



-1 q{q + a) pk \p — k\ + a 
For the calculation of the explicit CS term it is more convenient to use 



(A2) 

(A3) 
(A4) 



(A5) 
(A6) 
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' d. ^ = J- In (P + ^y + M' (A7) 

-1 q'^+p'^ + 2pqz + M2 2pq {p - qf + M"^ ' ^ ' 

1 pqz ^ ^p' + q^ + M^ {p + qf + M'' 

(iz— - — m — . (A8) 

-1 q^ +p^ + 2pqz + M2 Apq {p - q)^ + 



In pure QED this gives for the angular integration 



2a 1 2aMn(^ 



K(P,^)= I A^l .r. , — -3 li^l (A9) 



1 W (a + q^ q^ q 

2 a 1 ^ a + \k + p\ 

+ — In ■ ' ' 



max(p, A;) I A:^ — 7;^ I kp a+\k—p\ 



Q,2 in(i + «2 ^ ifc±ji^ 



+ , , . ."2 • (AlO) 



kp{—k + p)'^ kp{k+p)^ 



2. Explicit Chern— Simons Term 



The exphcit CS term, proportional to 



, fc^ dk /-i k-{k-p) 



k^ + MlJ-i {k-pf{\k-p\ + ay 

^'d^ , '^!;fVL . (All) 



can be calculated analytically. The angular integration gives 

Jo iq + ay\ Apq {p-q)^ + Ml)' ^ ' 

which we can calculate approximately by expanding the logarithm for p < g and p > q. 
Taking into account only the leading-order terms in min(p, q)/ max(p, q) gives 

3M|+p^ rP 2q'dq r-Q dq 2q^ 

^{M'i+p^Y Jo {q + af Jp {q + afMl + q^' ^ ' 

which can easily be calculated. Expanding the result for M± -C a gives 

3 Ml + p^ /p2 + 2ap 

3{Ml+p^y \ a+p 



+2 7 r ^ - - arctan + ^ In ^„ . A14 

However, we are only interested in the leading-order behavior in the infrared and ultraviolet 
region, so we finally get for p ^ a 

I{p) = 2 + 0{p/a) + 0(M±/a) , (A15) 

or p > a > M± 

I{p) = ^ + O(aVp') + 0{M^/p) . (A16) 
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APPENDIX B: UV EXPANSION OF THE HYPERGEOMETRIC FUNCTION 



1. Pure QED 

In pure QED, the ultraviolet boundary condition, Eq. |5^, leads to the requirement 

. 3 —a^ , 3 — 8 3 —a^ . , 

2Fi(a+, a_, -■ —^) + a 2^1(0+, a_, -; — ^) = 2^1(0+, a_, -; — ^) , (Bl) 



with a± = I ± ji^Nc/N — 1 and Ai'c = 128/(37r^). The derivative w.r.t. a of the hypergeo- 
metric function gives 

/ 3 — a^, 3 — a^, ^, 1 — a^. 

2Fi(a+, a_, -; — -) + a 2i^i(a+, a_, -; — r-) = 2i^i(a+, a_, -; — 5-) , (B2) 
2 m"' 2 ni'^ 2 m"' 

so we arrive at 

1 3 

2Fi(a+, a_, -; — -) = a+a_ 2^^1(0+, a_, -; — 5-) , (B3) 

2 m"^ 2 

where we have used 8/(37r^A^) = a+a.. To derive the behavior of the mass m close to the 
critical value, we expand Eq. (|B1| ) in a/m. Taking into account leading order only, we find 

r./ 1 r(l/2)r(a_ -a+) /«\-2a+ 

2 m"' 1 (a_)l (1/2 — a+) \m/ 

3 r(3/2)r(a_ -a+) /«\-2a+ 

.^iK, 2 ^ ^) = r(a_)r(3/2-a;) U + ^-^^ " ^^^^ 

We are interested in m/a, so we rewrite the boundary condition into 

(B6) 



m\'- _ r(l + ia;)r(a_)2(l - |a_) 



a 



r(l - iu;)r(a+)2(l - 



using the notation = y N^/N — 1. Since the absolute value of the RHS is equal to one, 
we can write it as eM<t>-niT)^ with 



r(l + >)r(a_)2(l-la_)l , (B7) 



= arg 

and n integer. This leads to the equation for m/a 

m -2mv + 20 

— = exp . (B8) 

a uj 

Next, we expand leading order in tu, using a Taylor expansion and the Euler constant 7 
r(l ± I c^) = 1 T icu7 + C(c^2) 

r(a±) = r(l/4) (^1 ^ (7 + Itt + 3 In 2) \uj^ + 0{u^) , (BIO) 

which leads to 
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l^ + hn2 + ]^]uj. (Bll) 



In this way we arrive finally at 

m 

— ~ exp 

a 



-2mx 1 2 
-7r + 3 ln2 + - 



2 7 



(B12) 



This expression is only valid for m ^ a, which means that n has to be positive. Using 
the effective potential, one can show that the largest value for mja gives the lowest energy, 
therefore the ground state corresponds to n = 1, and higher excited (oscillating) solutions 
are given by n = 2, 3, 4, . . .. 

For the chiral condensate we need to know the mass function at momentum p = a, m{a). 
Using the same expansions as above, we have 

^ 2^r(3/2)r(l - fg;) ( ay^'^^ , . \ .r.on 
m\a) = m —— — — - — - — + h.c. , (B13) 



^C(jr(a_)r(3/2 — a_|_) \m 
with m given by Eq. ( P8|) . Close to the critical coupling we can expand this in uo, using 



_ J(-7r+0) 

a 



^ -(cos0 + isin0) , (B14) 



which leads to 



m{a) „ r(l/2) fm\^l'^ ( f , vr 3 , \ ^ sin0\ , 
' ° I > I -H h - In 2 cos . (B15) 



a V{l/Af\aJ \\ ' 4 ' 2"" y " uj 

Finally, using the expansion for 0, we find 



a 7r(l/4)2Vay 
with m/a given by Eq. ( P12|) . 



m{a) ^ 64 1(1/2) fm^/^ ^ ^^^^^ 



2. Chern Simons Term: critical mass as function of 

The equation for as function of A^, Eq. |89|, is 

2Fi(a+, a_, -a^/M^) + a+a_ 2^^1(0+, a-, ^; -a^/M^) 

3 

-2a+a_ 2i^i(a+, O-, tj; -"V^c ) = • (B17) 

We can expand this as usual for small Mc 

/ a /r(-l/2)r(a_ - a+) r(l/2)r(a_ - a+) 

+ a+a 



\MJ \T{a^)T{-l/2- a+) ^ ~ r(a_)r(l/2 - a+) 

r(3/2)r(a_ - a 
r(a_)r(3/2 - a 



r(3/2)r(a_ - a+)\ , 
2a+a_ ^ ' . ±( +h.c. = 0, (B18) 
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which can be reduced to 



(1 + a+ + a+a_) + h.c. = , 



so we arrive at 



Mr 



a 



r(l + icj)r(a_)2(l + a_ + a+a. 
r(l - {uj)V{a^f{l + a+ + a+a_ 



Again, we can write this as 



Mr 



a 



exp 



-2n7r + 20 



UJ 



but now with 



arg[r(l + >)r(a_)2(l + a_ + a+a_)] . 



Finally we expand in a; 



'1 3, ^ 4\ 

~ i-vrH — ml cu, 

.4 2 21/ ' 



so we arrive at 



Mr 



a 



exp 



-2n7r 1 , 8 

+ -TT + 3 n 2 

UJ 2 21 



(B19) 
(B20) 

(B21) 
(B22) 
(B23) 
(B24) 



with n = 1,2,3, ..; the term —2nn/uj arises just as in pure QED, and also here there are 
infinitely many solutions for each value of N. The interesting one is the one corresponding 
to the ground state in pure QED, n = 1. 

3. Chern Simons Term: critical 6 as function of N 

Using the above approximation for Mc as function of A^, we can now calculate 6c as 
function of A^, see Eq. In general we have 

±TT^N 



16 



which we can expand 



M, (a+,a_, ^; -a^M^) - a+a_M, 2^1(0+, a_, | -a^M^) 



a 



A^r(l/2) fMA^/^ ( [MrY"^ r(a_ - 



"8(1 + 9^ 
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■ N' 



a 



a 



'1 - I a+) + h.c. 



We expand the F functions as before in uj and use 

'Mr\^'^ 



a 



(cos + « sm I 



which we have just derived, to get 



^ ^448 F(l/2) fMc\^'^ 



a 75 F(l/4)2 V a 
to leading order in using the expression for 0, Eq. P23 



(B25) 



(B26) 



(B27) 



(B28) 
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APPENDIX C: STATIONARY POINT OF THE EFFECTIVE POTENTIAL 



The effective potential given in section [V^ reproduces the SD equation for B± at the 
stationary point as follows. The variation with respect to B± of the effective potential is 
calculated as follows: 



Vo[B] 



6 



5B^{p)J (2^ 



In 1 + 



+ 2 



Blip) 



p^ + bUp) 



{2n 



and 



Vi[B] 



[p2 + 5^(p)]2 

2e2 [p^-Bl{p)] 



{2TTy [p^ + Blip)]^ J {27r)^ + Bl{k) 



X 



^2Driq) + ^]B^ik)T2^-^Do{q) 



Therefore we obtain 
6 



mip) 



V[B] 



2 [p' - Blip)] 
(27r)3[p2 + 5i(p)]2 

a{q) 



i^B^{p)-e'J 



d^k 



{27t)^ k^ + Bl{k) 



X 



2DT{q) + 



B±{k)T2^-^Do{q) 



(CI) 



(C2) 



(C3) 



This shows that the SD equation for B± is obtained from the stationary condition 
5V[B]/6B± = of the CJT effective potential V[B]. 
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FIGURES 



FIG. 1. Scaling of the fermion mass and the chiral condensates: the quantities lo ln(m(0)/a) 
and uj ln{{-ijjTp) / a'^) as functions of w = ^ N(./N — 1 for the hnearized equation (analytical solution, 
dashed line) and the full nonlinear integral equation (numerically, solid line) in pure QED. 

FIG. 2. The fermion mass function m{p) in pure QED: the analytical solution (dashed line) 
and numerical one (solid line) for N = 2 and = 3 where both the horizontal and the vertical 
axes are logarithmic. 

FIG. 3. The numerical solutions for wave- function renormalization functions, ^+(p) and 
A^{p), and mass functions, B^{p), B-{p), Be{p), Bo{p) aX 6 = 0.4 • 10~^ and m{p) at 6* = for 

= 3. (a) 1 — Aj^{p), 1 — A-{p), (b) mass functions in the IR region, (c) the same as (b) in the 
intermediate region. 

FIG. 4. The infrared values -B_|_(0) (dashed line) and -B-(O) (solid line) as functions of 
for some different values of 9, obtained numerically by solving the two sets of coupled integral 
equations. The dotted line indicates m(0) in the case of = 0. 

FIG. 5. The infrared values for i3+(0), i?_(0), -Be(O), -Bo(O), and the chiral condensate (V'V') 
functions of 6 for iV = 3, obtained numerically by solving the two sets of coupled integral equations. 
Note that the scale for the condensate is different. 

FIG. 6. The infrared values -B+(0) (dashed lines) and -B-(O) (solid lines) as obtained analyt- 
ically: (a) as function of for = 3, rescaled and compared with our numerical solutions, and 
(b) as function of A^ for some different values of 0, which should be compared with the numerical 
results plotted in Fig. ^ 

FIG. 7. Phase diagram in the (A^, ^)-plane: the critical line for the chiral phase transition 
obtained analytically and some numerical estimates for the critical point. 

FIG. 8. The chiral condensate and i?e(0) as functions of A^ for 9 = 10^^ and for pure QED 
(6* = 0), obtained numerically and analytically. Note that the scale for the condensate and -Be(O) 
is different. 

FIG. 9. The analytical solution of the linearized equation and the numerical solution of the 
full integral equation for A^ = 3 and 6 = 0.4 • 10~^. 
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